AD” AO 35  713  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CALIF  F/6  12/1 

ASYMPTOTIC  REPRESENTATION  OF  STIRLING  NUMBERS  OF  THE  SECOND  KIN— ETC (U) 
FEB  77  V E BLEICK.  P C WANG 

UNCLASSIFIED  NPS-53BL77021  NL 


END 

DATE 

FILMED 

3-7 


ADA035713 


NPS-53BL77021 

NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


Approved  for  public  release;  distribution  unlimited. 


Prepared  for:  j 

Office  of  Naval  Research  (Dr.  Bruce  McDonald) 

Statistics  and  Probability  Branch 
Arlington,  VA  22217 

COPY  AVAILABLE  TO  ODC  OSES  NOfl 
PERMIT  FULLY  LEGIBLE  PRODUCTION  i 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey*  California 


Rear  Admiral  Isham  Linder 
Superintendent 


J.  R.  Borstlng 
Provost 


ABSTRACT: 


The  distribution  of  the  Stirling  numbers  S(n,k)  of  the  second  kind 
with  respect  to  k has  been  shown  by  Harper  [Ann.  Math.  Statist.*  38 
(1967),  410-414]  to  be  asymptotically  normal  near  the  mode.  A new  single- 
term  asymptotic  representation  of  S(n*k),  more  effective  for  large  k*  is 
given  here.  It  is  based  on  Hermlte's  formula  for  a divided  difference 
and  the  use  of  sectional  areas  normal  to  the  body  diagonal  of  a unit 
hypercube  in  k-space.  A proof  is  given  that  the  distribution  of  these 
areas  is  asymptotically  normal.  A numerical  comparison  is  made  with  the 
Harper  representation  for  n-200. 


This  task  was  supported  by:  Contracts  No.  NR-042-286, 

NSWSES-56953, 

NISC-56969 


(/[>  z.fyLLi 


W.  E.  Bleick 

Professor  of  Mathematics  Emeritus 


Peter  C.  C.  Wang  ' 

Associate  Professor  of  Mathematics 


Approved  by: 


Released  by: 


jOfljgJtliaiiuisn  ^ 


IM 
1 IN 

■wiiortn 

1 “ 


Mathematics 


R.  R.  Possum 
Dean  of  Research 


NPS-53BL77021 


9 February  1977 


UNCLASSIFIED 


SECURITY  CLASS).  ICATION  OF  THIS  PAGE  fWAwi  Data  Knltnd) 


REPORT  DOCUMENTATION  PAGE 


4.  TITLE  (and  Submit) 


Asymptotic  Representation  of 
Stirling  Numbers  of  the  Second  Kindt 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


S.  RECIPIENT'S  CATALOG  NUMBER 


PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA  93940 


II.  CONTROLLING  OFFICE  NAME  ANO  AODRESS 

Office  of  Naval  Research  (Dr.  Bruce  McDonald) 
Statistics  and  Probability  Branch 
Arlington,  VA  22217 


NCY  NAME  t ADORES Mffaranl 


«.  DISTRIBUTION  STATEMENT  (ot  (ft/*  Ktporl) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  obetroct  on  toted  In  Block  20,  II  dltloront  from  Report) 


19.  KEY  WORDS  (Continue  on  roooroo  elde  II  nooooomr  **d  Identify  By  Block  » 

Asymptotic  representation 
Stirling  numbers  of  the  second  kind 
Bell  number 

Hermite's  formula  for  a divided  difference 


ABSTRACT  (Cbnllmf  on  ra*a?aa  m!d»  II  iMaaaaaqr  •**  Id—Itttr  br  Sfaa*  m 


....e  distribution  of  the  Stirling  numbers  S(n.k)  of  the  second  kind  with  respect 
to  k has  been  shown 'by  Harper  [Ann.  Math.  Statist.,  38  (1967),  410-414 ] to  be 


asymptotically  normal  near  the  mode.  A new  single-term  asymptotic  representa- 
tion of  S(n,k),  more  effective  for  large  k,  is  given  here.  It  is  baaed  on 
Hermite’s  formula  for  a divided  difference  and  the  use  of  sectional  areas 
normal  to  the  body  uiagonal  of  a unit  hypercube  in  k-space.  A proof  is  given 
that  the  distribution  of  these  areas  is  asymptotically  normal.  A numerical 
comparison  is  made  with 


DO  I i aZ'ti  1473  COITION  OF  I MOV  BS  IS  OSBOLCTS 

S/N  0 102*014"  4401  I 


- ' ' '4UVWW 


memm 


1.  Introduction. 

Previous  asymptotic  representations  of  Stirling  numbers  S(n,k)  of 
the  second  kind  have  been  of  two  types.  One  type  has  been  a complete 
infinite  series  expansion  as  given  by  Hsu  [1],  and  by  Bleick  and  Wang 

(2)  and  {3].  A second  type  has  been  the  single-term  representation  of 
S(n,k)  given  by  Harper  [4]  as  the  normal  distribution  approximation 

(1)  S(n,k)  ^ exp[-(k-u)2/2o2] 

o/*T 

where  the  mean  u 
numbers  B by 

(2) 

and 

(3) 

The  purpose  of  this  note  is  to  give  a new  single-term  asymptotic  re- 
presentation based  on  Hermlte's  formula  for  a divided  difference,  and  to 
compare  it  with  that  of  Harper. 

2.  Use  of  Hermite’s  formula. 

A Stirling  number  S(n,k)  of  the  second  kind  is  defined  as  the  kth 
difference  of  zn  at  z-0  divided  by  kl.  By  [5,p.l0]  we  find  that  this 
divided  difference  can  be  represented  by  a formula  of  Hermlte  as  the  re- 
peated definite  Integral 

1 t.  t,  , . 

(4)  S(n,k)  - / dt.  / dt,.../  K“A(du.n/du1k)dt. 

0 0 0 1 it 

where  u^«t^+t2+* We  imagine  that  t^»  t2t  t^  constitute  a set  of 


and  the  variance  o are  expressed  in  terms  of  the  Bell 


“ • Vi/B„  - 1 
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rectangular  Cartesian  coordlantes  and  Impose  an  orthogonal  transforma- 
tion of  coordinates  to  u^,  Uj»  • • > u^.  The  volume  of  the  space  over 
which  the  Integration  in  (4)  is  performed  Is  a portion  of  a unit  hyper- 
cube in  k-space.  If  we  allow  the  coordinate  u^  to  vary  along  the  body 
diagonal  of  the  hypercube  from  0 at  one  vertex  to  k at  the  opposite 
vertex,  the  sectional  areas  normal  to  the  diagonal  cut  by  the  hyper- 
plane Uj-tj+tj*. .+tfc  from  the  domain  of  integration  define  a positive 
function  gCu^k)  even  with  respect  to  the  argument  Uj-k/2.  We  take  the 
integral  of  g(u1,k)  to  be 

k 

(5>  / gCu^kJdUj^  - 1/k! 


to  agree  with  the  volume  of  the  space  over  which  the  integration  in  (1) 
is  performed.  We  drop  the  Uj^  subscript  henceforth.  Noting  that  g(u,k)-0 
for  k<u<0,  we  find  that 


{<3/ 2—  | u—3/2 1 ) 2/ 2 for  1/2  < |u-3/2|  < 3/2 
3/4  - (u-3/2;  for  1 <.  u <_  2 . 

Consideration  of  the  Laplace  transforms  of  (6),  (7)  and  (8)  suggests  that 


we  conjecture  the  Laplace  transform  of  k!g(u,k)  to  be 


for  all  k.  We  demonstrate  the  truth  of  this  conjecture  later 
ing  the  integration  in  (4)  over  the  variables  u„,  u„.  ...  u. 


Using  operation  82  of  [6,p.l0]  on  the  Laplace  transform  of 

u 

(11)  kl  j umg(u,k)du 

0 

we  find  the  mth  moment  of  the  k!g(u,k)  distribution  about  u-0  to  be 

(12)  lim  (-l)m(d/ds)B(l-e"8)k/sk  . 

s-K> 

It  is  now  easy  to  demonstrate  the  truth  of  the  conjecture  (9)  by  show- 
ing,  with  the  aid  of  the  multinomial  theorem,  that  (12)  is  the  same  as 


the  repeated  integral 


0 * 0 ‘ 

over  the  volume  of  the  hypercube. 


^odtl^dt2*  * ‘OW*  •+tk)m<ltk 


Use  of  (12)  and  (5)  shows  the  variance  of  the  k!g(u,k)  distribution 


to  be 


k/12  . 


Using  (14)  the  series 

as)  .,p(,,V/2) . i h-  ai pi  + + .. 

is  the  bilateral,  but  not  s multiplied,  Laplace  transform  of  the 


normal  distribution 


(l/a/2iT)exp(-t2/2o2) 


according  to  [7,p.2|].  The  corresponding  series  for  (9)  multiplied  by 
|f0  / 2 

e , or  the  bilateral  Laplace  transform  of  k!g(u,k)  shifted  left  by 


k/2,  is 
(17) 


2 /,  v 2 


(2/s)ksinhks/2  - [1  + + ..]k  . 


The  dominant  k power  term  in  the  coefficient  of  (s2/4)n  in  (15)  is 
kn/6nnl,  and  may  be  shown  to  be  the  same  in  the  expansion  of  (17) 
by  the  use  of  the  recurrence  formula  6.361  of  [8,p.ll9].  This  proves 
that  the  klg(u,k)  distribution  is  asymptotically  normal  as  le**».  It  is 
rsmarkable  that  the  normal  distribution  should  arise  in  the  purely 
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geometrical  context  of  sectional  areas  normal  to  the  body  diagonal  of 

a hypercube  of  high  dimension. 

% 

On  replacing  klg(u,k)  in  (10)  by  its  Gaussian  normal  approximation 
of  mean  y-k/2  and  variance  o2-k/12  we  find 

(18>  s<n»k>  * -7rfi>  I un‘kexp[-(u-k/2)2/2o2]du 

/3  k 2 

(k/2-at)n-'ce-t  /2dt  . 

*■  -CO 

3.  Numerical  example. 

Table  1 compares  the  exact  values  of  S(200,k)  with  the  asymptotic 
approximations  computed  from  the  single-term  representations  (1)  and 
(18).  Harper's  representation  (1),  which  uses  B^-. 62475  10276’ 
y-49.975  and  0-3.0551,  gives  an  excellent  fit  near  the  mode  (k-50) , 
but  (18)  gives  a much  better  fit  for  large  values  of  k. 

Table  1.  Values  of  S(200,k) 


k 

Asymptotic 
from  (1) 

Exact 

Asymptotic 
from  (18) 

2 

.23135  10222 

.80347  106 

.69244  10126 

40 

.39504  10273 

.24458  10273 

.42658  10273 

50 

.81579  10275 

.81493  10275 

.15285  10277 

60 

.37452  10273 

.53533  10273 

.29658  10274 

100 

217 

.49065  104  ' 

.22839  10235 

235 

.27994  10 

150 

.13938  1043 

.30251  10143 

.30441  10143 

199 

.16955  10“241 

.19900  105 

.19900  105 
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